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Abstract 



The Direct and the Inverse Scattering Problems for the heat operator with a po- 
tential being a perturbation of an arbitrary N soliton potential are formulated. We 
introduce Jost solutions and spectral data and present their properties. Then, giving 
the time evolution of the spectral data, the initial value problem of the Kadomtsev- 
CNj , Petviashvili II equation for a solution describing N solitons perturbed by a generic 

smooth fast decaying potential is linearized. 

m 

. 

• 1 Introduction 

. The Kadomtsev-Petviashvili equation in its version called KPII 

(u t - 6uu Xl + u XlXlXl ) Xl = -3u X2X2 , (1.1) 



where u — u(x,t), x — (xi,X2) and subscripts xi, X2 and t denote partial derivatives, is 
a (2+l)-dimensional generalization of the celebrated Korteweg-de Vries (KdV) equation. 
The KPII equations, originally derived as a model for small-amplitude, long- wavelength, 
weakly two-dimensional waves in a weakly dispersive medium [1], have been known to be 
integrable since the beginning of the 1970s [2,3], and can be considered as a prototypical 
(2+l)-dimensional integrable equation. Its integrability results from the existence of the 
Lax pair 

C(x,d x ) = -d X t+% 1 -u{x), (1.2) 

xi 

A{x,d x ) =49^ - 6ud Xl -3u Xl -3 / dx[ u X3 (x[, x 2 ), (1.3) 



so that (1.1) is equivalent to the compatibility condition 

C t = [C,A]. (1.4) 

Operator (1.2) defines the well known equation of heat conduction, or heat equation for 
short. The corresponding Inverse Scattering Theory (1ST) was developed in [4-8] in the 
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case of a real potential u(x) rapidly decaying at spatial infinity. This theory is based on the 
integral equation 

$(x,k)= e - 4k ^- k2 ^+ J d 2 ygo(x-y,k)u(y)<S>(y,k), (1.5) 

for the Jost solution. Here k e C is the complex spectral parameter and Q n (x 7 k) is the 
Green's function of the bare heat operator — d X2 + d Xl , given by 

g Q (x,k) = -^l^ie^k^-k 2 ^ J ds e( s { s + 2k u )x 2 ))e- ts{x ^- {s+2 ^\ (1.6) 

where 9 is the Heaviside function and integration is along the whole axis. However, class of 
such potentials is not the most interesting one, since the KPII equation was just proposed 
in [1] in order to deal with two dimensional weak transverse perturbation of the soliton 
solution of the KdV. This solution, as well as multisoliton ones, is known to have non 
decaying, ray asymptotic behavior at space infinity (see, e.g., [9]- [14]). In order to illustrate 
this behavior, let us notice that equation (1.1) is invariant under the Galileo transformation, 
i.e., if u(xi,X2,t) is a solution, then 

u(xi, X2, t) — u(xi + [1X2 — 3[i 2 t, x 2 — Q[it, t), (1.7) 

obeys this equation as well for any real constant [i. Then the ray asymptotic behavior 
follows because any solution u(xi,t) of the KdV equation solves KPII. This behavior, [9], 
causes the main difficulty for construction of the spectral theory of the operator (1.2) with 
a potential 

u(x) = un(x) + u'(x), X = (xi,X2), (1.8) 

where un(x) is some N soliton potential and u'(x) is a smooth and rapidly enough decaying 
function of its variables, that can be considered a perturbation of the soliton potential. 
Indeed, let us consider the case of the pure one soliton potential, u(x) = u\(x\). Then the 
Jost solution has the form $(x, k) = e jk:E i+ k X2 x(xi,k), where x( x i,k) is a function of 
only one spatial variable. Inserting this function into (1.5) we see that the integral with 
respect to 7/2, as follows from (1.6), is divergent. Therefore, already for the simplest case of 
the potential of the type (1.8) the standard integral equation (1.5) for Jost solutions is ill 
defined, as its kernel does not exist. 

Our aim here is the construction of a spectral theory of the KPII equation that also 
includes solitons, as already successfully done for the KPI equation [15]. Following the 
framework of that article we develop a generalization of the standard 1ST: the scattering on 
a nontrivial background. More precisely, we take, from the literature, the known generic ex- 
pressions of the pure N soliton potential Un(x) and the corresponding Jost solution $at(x, k) 
of the operator Cn, i.e., C in (1.2) with u = ujy. Such potentials and Jost solutions can be 
constructed either by means of rational similarity transformations of the spectral data, or 
by the r- function approach, or by twisting transformations (sec [9,11,16] and Sec. 2 here), 
as all these procedures lead to the same result, as it was shown in [17]. Then we define 
the Jost solution of the operator C with potential u (1.8) by means of the following integral 
equation 

$(x, k) = $ N (x, k) + J dy Q N (x, y, k)u'(y)$(y, k), (1.9) 

where Q n{x, y, k) is the Green's function of the operator Cn, i.e., solution of the differential 
equation 

(-d X2 + d 2 Xl - u N (x)) G N {x,x', k) = S(x - x'), (1.10) 
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and its dual, fixed by the condition that 

G N (x,x' ,k) = e *k(*i-*i)+k a (* a -*i) g N ( x ,x' ,k) (1.11) 

is bounded with respect to the variables x, x' € K 2 and k € C and has finite limits at infinity 
( "boundedness condition"). In contrast to (1.5) the integral equation (1.9) is well defined 
thanks to this condition either under small norm assumptions on u'(x), or under assumption 
of decaying of this perturbation faster than some exponent linear in x. In addition, it is 
natural to expect that the proof of unique solvability of (1.9) could be performed in analogy 
to the case of a decaying potential, see [6-8]. 

Here, we assume this unique solvability and study the properties of the total Green's 
function, Jost solutions and scattering data of the perturbed potential u(x), exploiting the 
fact that their properties are inherited, via equations of the kind (1.9), from properties of 
these objects in the pure N soliton case. This approach was already developed in [18] for 
the perturbation of the one soliton potential (see also [19], where some class of potentials 
non decaying in one space direction was considered), while the generic case of the iV-soliton 
potential was left open just because the Green's function obeying boundedness condition 
was unknown. This problem was solved in our article [21] (preliminary version was given 
in [20] ) , where by means of the extended resolvent approach we derived the Green's function 
Qn{x, y, k) and presented its properties in detail. In particular, we showed that this function 
is discontinuous at some points on the k-plane, determined by the soliton parameters, and 
we derived a set of auxiliary Green's functions that enable the control of the discontinuities 
of Qn{x, y,k) and Jost solutions, see Sec. 3 below. 

The article is organized as follows. In Sees. 2 and 3 we give necessary details for the 
pure N soliton potential, Jost solutions and Green's function. The Direct problem for 
the perturbed potential (1.8) is given in Sec. 4. In Sec. 5 we investigate properties of 
the corresponding Jost solution and scattering data and derive the Inverse problem. Time 
evolution of the scattering data is given in Sec. 6. Concluding remarks and some future 
developments are presented in Sec. 7. 



2 Heat operator with multisoliton potential and its Jost 
solutions 

Soliton potentials (see [16] and [9-14,17,22] for details) are labeled by two numbers (topo- 
logical charges) N a and Nb, which obey condition 

N a ,N b >l. (2.1) 

Let 

Af = N a + N b , so that M > 2. (2.2) 
We introduce the M real parameters 

Ki < k 2 < ■ ■ ■ < km, (2.3) 

and functions 

K n (x) = K n x x + nlx 2 , n=l,...,N. (2.4) 

Let 

e K ^ = diag{e K "^}^ =1 (2.5) 
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be a diagonal Af xAf matrix, let T> be a AfxNi, real constant matrix, with at least two 
nonzero maximal minors, and let P' be a constant N a x Af matrix that, like the matrix 
T>, has at least two nonzero maximal minors and that is orthogonal to the matrix T> in the 
sense that 

£>'£> = 0, (2.6) 

where the zero in the r.h.s. is a N a x Nb matrix. Let us also introduce the two "incomplete 
Vandermonde matrices" (Nb x Af and Af x N a , correspondingly) 



V = 



/ 1 

Kl 



1 



AW 



( 1 ... 



N a -1 



V 



(2.7) 



Vl ... nfr-i / 



In analogy to (2.6) these matrices obey the orthogonality condition 

V 7 V' =0, 

where 7 is the Af x Af constant, diagonal, real matrix 



M 

7 = diag{ 7 „}^ =1 , 7„ = Yl («« - Kn') _1 - 

n'=l 
n T^n 

Then, the soliton potential is given by any of two representations 
u N (x) = -2d 2 Xi \ogr N {x) = -2d 2 Xi \ogr' N {x), 
where the r-functions are expressed as 

t n (x) =det(Ve K ^V), T , N (x)=det(v'e- K ^ 1 V') . 



(2.8) 



(2.9) 



(2-10) 



(2.11) 



For the Jost and dual Jost solutions (solutions, respectively, of the heat operator (1.2) and 
its dual) we have 



-ikii-k x 2 



TtJxM ^ ( k)=e ikx 1+ k'x a Wg;jf) (212) 

T N (X) T N {X) 



where r $jv (2, k) and (x, k) are given in terms of the T-functions by means of the Miwa 
shift: 



T<s> N (x,k) = det(Ve K{x) {K + ik)V), T* N (x,k) = det ( V 

and where we introduced one more Af x Af diagonal matrix 

k + i k = diag{«„ + i k}^. 



,K{x) 

K + ik 



V 



(2.13) 



(2.14) 



Analogous relations are valid in terms of rjy. Thanks to these definitions we have the 
following conjugation property of the Jost solutions: 



$ N (x,k) = $ N (x,-k), ^ N (x,k) = * N (x,-k). 



(2.15) 
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In order to study the properties of the potential and the Jost solutions, it is convenient 
to use the representations for the r-functions, which follow from Binet-Cauchy formula for 
the determinant of a product of matrices, i.g. for the r-function we have 



T N {x) 

t*„ (x, k) 



— ]T V({n i })V({n i })He K 



.,(*) 



{n,}=l 
Af 



1=1 
N b 



1 _ _ _ _ 

— ]T V({n l })V({n l })l[( Kni + l k)e K 



{",}=1 
N 



1=1 

2± e ^n ; (x) 



{n,}=i 



1=1 



+ ik' 



(2.16) 
(2.17) 
(2.18) 



where we use notation 
1 



v(K» - 



K 



N b -1 



K 



N b -1 



z>({M) = 



^W, ,1 



V 



ni,N b 



V 



(2.19) 



7ljV b ,N b 

,riN b } stands for not 



for the maximal minors of matrices V and V and where {rii} = {rii 
necessary ordered sets of Nt indexes from the interval l,...,Af. 

We recall that the maximal minors of a matrix satisfy the Pliicker relation, i.e., for any 
subsets {mi} and {r^} of indexes running from 1 to N and arbitrary j <G {1, . . . , N^} 



V({m i })V({n i }) = 

N„ 

= ^ £ ) ( m i; ■ • • , m s -i,nj,m s+ i, m Nb )V{ni,. . . , n,_i, m s , n j+1 , n Nb ) 



(2.20) 



The matrices V and D are not in one-to-one correspondence with the potential u(x). Indeed, 
from (2.10) and (2.11), one gets that the potential is invariant under the substitutions 



(2.21) 



where v and v' are respectively arbitrary real constant, nonsingular Nb x iVj, and N a x _/V a 
matrices. Thus under condition (2.3), the (N a , Ni,)-so\iton potential is parameterized by a 
point on the Grassmanian Gr^w/ if the first representation in (2.11) is used, or by a point 
on its dual parametrization (i.e., the Grassmanian GiN a ,M) ^ the second representation in 
(2.11) is used. 

Regularity of the potential u(x) on the x-plane is equivalent to the absence of zeros of 
t(x). It is clear that it is enough to impose the condition that the matrix T> is Totally Non 
Negative (TNN), i.e., that 



T>(m, . . . ,riN b ) > 0, for all 1 < m < . . . < riN b < N ■ 



(2.22) 



It is known that in the classification of the TNN matrices a special role is played by the 
solid minors [23] , i.e., minors of the kind T>(n, . . . , n + Nb — 1) for all n = 1, . . . , Af where we 
assume that the indexes are defined mod A/", so that thanks to (2.2), say, n + Nb = n — N a 
for n > N a . In [10] we proved that if the coefficients 



V{n, ...,n + N b -l) V(n, ...,n + N b -l), 



(2.23) 
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are strictly positive for all n = 1, . . . , A/" (again with indexes defined mod N) the function 
t{x) has the following asymptotic behavior 

/n+N b -l \ 

x^Aoo: T(x) = (z n + z n+1 e KN >+^- K ^+o(l))cxp\^ ]T K^J , (2.24) 

/n+N b -l \ 

x^oo: t(x) = (z n + o(l))cxpi K i( x )Y ( 2 - 25 ) 

Here r„ denotes the asymptotic rays on the x-planc 

f x 1 + ( Kn + K +Nb )x 2 bounded = 

1 (Kn+W b - k™)z 2 -> -oo. v 7 

so that there are N a rays in the direction x 2 — > — oo and iV& rays in the direction X2 — > +oo. 
As well, cr„ denotes the sector swept out by rotating anticlockwise the ray r n up to the ray 
r n+ \. These sectors are nonintersecting and cover the whole x-plane with the exception of 
rays. 

It is clear that the condition 

z n > (2.27) 

is sufficient for having nonsingular asymptotics of the potential. In the case of a TNN matrix, 
this condition is equivalent [24] to the condition that all maximal minors of the matrix V 
are positive, i.e., I? is a totally positive (TP) matrix. 

A special role in the construction below is played by the values $Ar(x, in n ) of the <&n(x, k) 
at points in n and the residuals of \E , Ar(x,k) at the same points. In terms of these discrete 
values we define a A/"- row 

®n, k ( x ) = {$n,i{x), $iv,jv(a;)}, ®N,n( x ) = ®n{x, iK n ), (2.28) 
and a A"-column 

^nA x ) = {^n.i(x),---,^n,m(x)} t , ^nA x ) = ~i res *(z,k), (2.29) 

k— in n 

n = 1, . . . ,Af. In [22] (in a bit different notation) we proved that 

$ N , K (x)V = 0, (2.30) 
T>'* NtK (x) =0, (2.31) 

and that the Jost solutions obey the Hirota bilinear identity 

N 

J2®nAx)VnAx')=0- (2-32) 

n=l 

Thanks to (2.15) we have that 



$N,n( x ) = ®N,n(x), ^nA x ') = $nA x '), ( 2 - 33 ) 

for all n. 

In [9,10,17,22] detailed properties of the Jost solutions and their discrete values are given. 
Here we mention only that the functions e lkxi+k X2 <&n(x, k) and e -* ka; i- k X2 <$) N (x, k) 
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have bounded asymptotics on the x-planes, function e lka:i+k X2 <$?n(x, k) is a polynomial 
with respect to k of order Nf,, and function e - lka; i-k X2 ty N (x, k) is a meromorphic function 
of k with simple poles at all k = in n and of order — Nf, at infinity. More exactly, we have 
the following normalization 

lim {ik)- N »e ikxi+]52x2 $ N (x,k) = 1, lim {ik) Ni e- ikxi -^ x '^ N {x,k) = 1, (2.34) 

k— foo k— >oo 

so that the potential can be reconstructed as 

u N (x) = -2 lim (ik)- N > +1 d Xl (e i * Xl+ * X2 & N (x,k)) = 

k— >oo 

= 2 lim (ik) JV » +1 a xi (e- ika,1 - k2xa * JV (a:,k)). (2.35) 

k->oo 

Analyticity properties of the Jost solutions with respect to k impose conditions on their 
discrete values. Indeed, taking notations (2.4), (2.5), and (2.9) into account, we can write 
that 

e ik, 1+ k» „g ( k) = f[ {Km + ik) j2 -rne-^ w^„(») (2 36) 

On the other side we know that the l.h.s. is a polynomial of the order Nf,, normalized 
according to (2.34). Thus, since all higher powers in the above equality must cancel out, we 
obtain that 

<S>nAx)i^ K[x) V - { 0,...,0 , {-l) N ") . (2.37) 

An analogous formula for &n,k(x) can be derived directly, in a similar way, or by recalling 
(see [22]) that 

<M*,k)= (2.38) 
ft K + ik) 

n=l 

where square brackets in the r.h.s. means that one has to perform in $jv(— x,k) the substi- 
tutions V — > 7l?' T , N a <H> JVj,, while the k„'s remain unchanged. Precisely, we have from 
(2.38) that *„(x) = {-1) M ^ n [^ Ntn {-x)\, and, then, thanks to (2.37) we get 

V e K ^ N , K (x) = (0, . . . , (-1)^) T . (2.39) 

Equalities (2.30) and (2.37) define the discrete values &N,n(x) and then the Jost solution 
$jv(x,k) itself by (2.36). The same is valid, of course, for the dual Jost solution ^jv(x,k), 
where the defining equations are (2.31) and (2.39). The potential ujv(x) can be reconstructed 
by (2.35). 

In order to simplify relations below, we introduce operatorial notation of the kind 

l N $ N (k) = 0, *iv(k)1c N = 0, (2.40) 
for the heat equation and its dual and 

l N g N (k) = g N (k)t N = i, (2.4i) 

for the equation (1.10) and its dual. Thus ~L jv denotes the operator Cn applied to the 
x variable of the kernel Qn{x, x',k) and Cn denotes the operator dual to Cn (i.e., £%) 
applied to the x' variable of the same kernel. 
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3 Green's functions of the pure soliton potential 

3.1 Green's function g N (x, x') 

In [21] by means of the extended resolvent approach we derived the Green's function of the 
heat operator with generic TV-soliton potential obeying condition (2.27) (for a preliminary 
version see [20]). In terms of the Jost solutions it can be written in the form 

g N (x 7 x\k) = - sgn ^ x ' 2) Jd s e(( s 2 -kl)(x 2 -x' 2 ))<s> N (x, s + ik^ N (x', s + ikz)- 



- 6{x' 2 - x 2 ) 0(ka -K n )^ N!n {x)^ Nin (x'). 



(3.1) 



We proved there that function Gn(x, x' k') defined in (1.11) is bounded with respect to the 
variables x, x' € K 2 and keC and has finite limits at infinity. We also derived the auxiliary 
Green's function 



9n(x,x') 



0(x 2 



2tt 



— |y dae~ ia ( Xl ~ x '^~ a2 ( X2 ~ x ' 2 ^ + 



TV . , 

+ lJ2^ N . n (x)^ N , n (x') / — e -i<x(xi-x'i+2Kn(x 2 -x' 2 ))-a 2 (x 2 -x' 2 ) I ( 3 2 ) 

n=l ) 

where integrals are understood in the sense of the principal value (in [21] this function was 
denoted as 9 + ). Function Gn(x, x' , k) is continuous with respect to the variable k <G C with 
exception of the points k = im, . . . Discontinuities at this points can be described 

by means of the auxiliary Green's function (see below). This function can be calculated 
explicitly as 



/ r\ @{x 2 X 2 ) 

g N (x,x ) = , , ==exp 



[x x -x\f 



2 v /ir(x 2 - x' 2 ) 



A(x 2 - x' 2 ) / 

X ^ X ~\~ QitCfi^X^ 



erf 



4) 







(3.3) 



where the error function is defined by 

erf(x) 

so that 



2 

V5F 



dye y , 



lim erf(x) = 1, lim [erf(x) — l]xe x = —=. 



(3.4) 



(3.5) 



Taking into account (2.32), we subtracted 1 in the square brackets of the r.h.s. for making 
explicit the good behavior at large x of the last term in the r.h.s. thanks to the first limit 
in (3.5). To specify this behavior we consider 



6>(x 2 - x' 2 ) 
2 a /tt(x 2 - x' 2 ) 



exp 



(x x — x\ + 2 k 3 (x 2 — x 2 )Y 
4(x 2 - x' 2 ) 
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- ^ 2 ^ jr e^"( a; )$ A r^(x)e- K "( a; '^ 7 v,„(x , )e( K "- k3 )( a: i- a; 'i+ 2K "^-^ 



2JJ X 



x e 



n=l 

-(K n -k a ) 2 (x 2 -x' 2 ) 



erf 



2K n (x2 — X2) 



(3.6) 



When x — > 00 the first term in the r.h.s. is decaying in all x directions and for any 
at least as (x 2 — x^ 1 / 2 . Let us consider the second term taking into account that all 
e Kn ^ x ^N,n{x) and e~ Kn ^ x ^ N,n{x') are bounded for all x and have finite limits at infinity. 
The first exponent under summation sign can grow linearly with respect to the combination 
x x — x[ + 2n n {x 2 — x' 2 ). But for generic situation this is depressed by the behavior of the 
second exponent and erf in the second equality of (3.5). Only the special case, in which the 
asymptotic behavior of x is such that for some n the argument of erf is bounded, is left. 
But then the first exponential factor is depressed by the second one if ks ^ k„. If kcj = n n , 
then both these exponents are absent and again we have boundedness of the r.h.s. in (3.6). 
This proves that for any kg g 1 the product 



e -k 3 K-<)-k|(z 2 -4) g ^ (a ,y ); 
is bounded. Let us mention also that the Green's function g^ is real, 



g N (x,x') = g N (x,x'), 
as follows from (2.15), (2.33), and (3.4). 

3.2 Green's function Q N (x,x' 
By definition function Qn{x,x' ,k) obeys 



(3.7) 



(3.8) 



g N (x,x',k) = g N (x,x',-k) = g N (x,x',k), 

thanks to (2.15), (2.33), and in terms of distributions we have by (3.1) 



dg N (k) _ sgnkgj 



<9k 



2tt 



$jv(-k) ® *jv(-k), 



(3.9) 



(3.10) 



where for shortness we omit space variables and use the notation of the direct product in 
the standard way: ($ N (-k) (g) ^ N (-k))(x,x') = $ N (x,-k)fy N (x' ,-k). In [21] we proved 
that the Green's function (3.1) can be written in the form 



g N (x,x',k) = g N (x,x')+ 



1 



IM 



— / da 

27r ^-|k !R | 

N 



$n(x, a + i k 9 )*Ar(x', a + i k a ) - ^ 



kg —n m - ia 



+ ^ a(k-iK m )^N,m(x)^ N,m(x'), 



m—1 



where wc introduced the function 



«(k) = iarccot^ 

71" I k R I 



+ 



(3.11) 



(3.12) 
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Function a(k) satisfies the following relations 



a(k) = a(k) = a(-k), < a(k) < 1, (3.13) 
a(ie ia | k |) = l^i, -7T < a < tt, (3.14) 

Since the integrand in (3.11) has no singularities, we see that, indeed, the Green's function 
Gn(x, z',k) is discontinuous at all points k = ik„ and has finite limits at these points 
depending on the limiting procedure (see (3.14)). In order to control these discontinuities it 
is convenient to introduce an auxiliary Green's function that in the case of the heat equation 
(see [21]) coincides with the value of Gn(x,x' ,\c) at kg = 0. Since the term in the second 
line of (3.11) tends to zero when kjj — > we have by (3.11): 

Gn{x, x', i kcj) = g N (x, x') + ^ ^N, m (x)^N, m (x')9(K m - k 9 ), (3.15) 

m— 1 

which is defined for any value of kg ^ n n , n = 1, . . . , J\f. It is discontinuous at all points 
kcj = K\, . . . , njsf and piecewise constant function of kcj otherwise. In the case where kcj < 
Ki, or kg > k^/ we have thanks to (2.32) that Gn(x, x', i kcj) = g^{x,x') and 

lim Gn(i^) = lim Gw(i k a ) + <$> N , n <g> ^ N ,n- (3.16) 
Because of (2.33) and (3.8) this function is real 



G N (ik*&)=G N (iksi), (3.17) 
and taking its discontinuity into account we denote 

G Nn = lim lim &v(k), n = l,...,M. (3.18) 

k Q ^K„+0ksH^0 

Then by (3.15) we have that 

M n 

G N ,n =9N+ ^N,m{x)^N,m(x') = 9N~Y1 ^N,m{x)^N,m{x'), (3.19) 

m— n+1 m— 1 

where indexes are defined mod A/", as it was mentioned after (2.22), and where the second 
equality follows by (2.32). Thus 

G N ,n-l = G N ,n +®N,n ® *JV,n. (3-20) 

and then for any k in a neighborhood of a point in n and any n — 1, . . . , N we have: 

5jv( k ) = £jV,n+ a ( k - iK n) $ JV,n ® *JV,n + o(l), k = IK„. (3.21) 

4 Perturbation of iV-soliton potential 
4.1 Main definitions 

Let us now consider the heat operator (1.2) with an TV-soliton potential perturbed by adding 
to it a real function u'(x), smooth and rapidly decaying at space infinity, see (1.8), i.e., 

C(x, d x ) = C N (x, d x ) - u'(x). (4.1) 
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Below we use the operatorial notation C = C N —U', where U' denotes the multiplication 
operator with kernel 

U'(x,x') = u'{x)5{x-x'). (4.2) 

Thus, we write the integral equation (1.9) for the Jost solution of the heat operator C in 
(4.1), and for its dual, as follows 

$(k) = <Mk) + g N (k)u'Mk), *(k) = <Mk) + *(k) V g N (k). (4.3) 

In fact we introduce and study here a more general object: the total Green's function 
G(x, x' , k) defined for any value of the complex spectral parameter k by means of one of the 
following integral equations: 

0(k) = g N (k) + g N (k)u' e?(k), g(k) = g N {k) + g(k)u' g N {k). (4.4) 

As it was mentioned in the Introduction, we assume here that these equations are 
uniquely solvable and give the same solution at least for smooth u'(x) obeying some small 
norm conditions, or decaying fast enough at space infinity, e.g., faster than some decaying 
exponent linear in x. We also assume that this function inherits properties of the pure 
solitonic Green's function <?Ar(k) by means of (4.4). In particular, g(x,x',k) obeys the 
differential equations 

tg{k) = g{k)t = i (4.5) 

and, thanks to (2.41), boundedness property (1.11). This function is continuous, continu- 
ously differcntiable with respect to the complex spectral parameter k, for all values of this 
parameter different from im, . . . , in^f- Because of (3.9) and reality of the potential v! we 
have the conjugation property 



g( x , x',k) = g(x, x', -k) = g(x, x ', k). (4.6) 

In terms of the Green's function we define the Jost and dual Jost solutions as 

Hk) = g(k)t N ^ N (k), *(k) = y N (k)t N g(k), (4.7) 

where, say, the first equation is a short form for 

$(x,k) = J dx'{£%(x',d x ,)g{x,x',k)}<t> N (x',k). 

Then, thanks to (4.1), (2.40), and (2.41), the differential equations 

Z$(k) = 0, *(k)£ = 0, (4.8) 

hold and, therefore, we indeed can consider $(x, k) and ^(x, k) as the generalization of the 

Jost solutions to the case where the perturbation u'(x) is different from zero. Because of 
(4.1) we can write (4.7) also in the form 

$(x, k) = $ N (x, k)+ J dx' g(x, x', k)u'(x')<S> N (x', k), (4.9) 

*(x, k) = * N (x, k) + J dx' g(x', x, k)u'(x')y N (x', k), (4.10) 
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Thanks to the above assumptions both these solutions are continuous, continuously differen- 
tiable function of k e C for all values of this parameter different from in\, . . . , i/t/v- Because 
of the asymptotic conditions (1.11) and (2.34) we have that 

lim (ik)- JV *e <kxi+k2xa $(a;,k) = 1, lim (ik) N "e~ lkxi - ]i2 X2 ^(x,k) = 1, (4.11) 

k— >oo k— >oo 

and thanks to the reality of the perturbation u'(x) and (4.6) we have the conjugation prop- 
erties 

$(x,k) = $(x,-k), *(avkj = *(x,-k). (4.12) 

4.2 Behavior of the Green's function at points of discontinuity 

The Green's function Qjy(x,x' ,k) has a well defined limit at ksR = 0. Correspondingly, 
the limit Q(ik<^) at ksR = of the total Green's function obeys the integral equations 
C/(ikcj) = g N (ik^) + ^Ar(ikcj)[/'5(ik Q ) and £/(ik Q ) = £jv(ikc>) + Q{i k^)U' Q N (i k 9 ), 
obeys the differential equations ~E Q(i kcj) = C/(i kcj) C = I, and is a real, piecewise constant 
function of kg. Like (?_/v(*kcj) it has finite limits from the right and from the left at points 
k$j = K n . In analogy to (3.18) we denote 

g n = lim 0(ikc»), n = l,...,AT, (4.13) 

k 3 — 

that by (4.3) obeys the integral equations 

G n = G N , n + G N , n U'Q n , G n = G N , n + G n U'G N , n - (4.14) 

These values of the Green's function can be related by means of the "dressed" version of 
equation (3.20). Say, by the first equation above we get 

Q n -\-Q n = [GN,n-i-s N , n ]{i + u'g n ) + g Ntn _ 1 u'[g n _ 1 -g n \. 

Thanks to (4.1) and (4.14) we can insert here I + U' g n = ~£nGti- Then using (3.20) for 
the difference Gn n-i~ Gn n we derive 

Gn-i -G n = *jv,« ® *n,u~Z n g n + g N ^ u'\g n _i - g n \. 

We see that the difference g n -i — Gn obeys an integral equation with the same kernel as in 
(4.14) (for n — > n — 1) but with a different inhomogeneous term. Thanks to the equality 

g N,n-i C-n = I this term can be obtained from g N n _ 1 by applying operation C N$N,n®^ n 
to (4.14) from the right. Thus under condition of unique solvability of this equation we arrive 

to 

g n _ 1 -g n = g n _ 1 £ N $N, n ®* n , (4.15) 

where we introduced the auxiliary solutions of the heat operator and its dual (cf. (3.18)): 

$ n = g n t N$N,n, (4.16) 

* n = V Nin l N G n , (4.17) 

n € [1,A/], which are real and satisfy = = $„ C . They can be considered as the 

discrete values of the Jost solutions, if we notice that thanks to (4.7) and (4.13) we have the 
following generalizations of (2.28) and (2.29): 

$„(x)= lim $(x,ik Q ), (4.18) 

k 3 — >re„+0 
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*„(x) = lim (k 9 -K n )*(a;,ikc f ). (4.19) 

kg— >K n +0 

The first factor in the r.h.s. of (4.15) is not of the kind (4.16), so we apply operation 
£ N$N,n to (4.15) from the right, that gives 

(g n _ 1 t N * N , n ){x) = ^-, (4.20) 

where we introduced the constants 

c„ n=l,...,7V. (4.21) 

Thanks to (4.1), (4.8), and (4.17) it can be written in the following different forms 

c n - {^n^n Q n = (? Nin ? N * n ) = = (* n U'$ Nin ) , (4.22) 

where, say, the last equality is a short form for 

c n = j dxu'(x)y n (x)<S> N Jx). (4.23) 

Because of the boundedness property (4.11) and thanks to (2.4) and (4.18) we see that in 
analogy with the unperturbed 7V-soliton case products 

e- Kn{x) $ n {x) and e K " (x) ^ n (x) are bounded for all x (4.24) 

and have finite limits at ^-infinity. Thus, the integral in (4.23) exists under some small 
norm condition on u'(x), as discussed above. Then since both sides of (4.20) exist and 
$„(x) cannot be identical zero, we get that c„ ^ 1, and thus 

c n < 1, (4.25) 

because in the limit of the zero perturbation, u'(x) — > 0, all these constants tend to zero. It 
is also clear that all these constants are real. 

Finally, by (4.15) and (4.20) we get the following generalization of (3.20) for the perturbed 
case: 

y n -\=Qn-\ : , (4.26) 

1 - Cn, 

so that, summing over n (under mod AA-condition), we get 

E%^=0, (4.27) 

n=l 1 C ™ 

i.e., a generalization of (2.32). 

On the other side, summation of (4.26) for some interval of indexes in [1,7V] gives the 
equality 

g m = g n + J2 (4 - 28) 

l=m+l 1 

that under mod TV-condition is valid for any m, n e [1,7V]. 
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In analogy to the above we can consider the difference <?(k) — Q n . Thanks to the first 
equalities in (4.4) and (4.14) we have 

S(k) - g n = [g N (k) - g N jl N g(k) + g N , n u'[g(k) - g n ], 

so that for k in a vicinity of some in n , n= l,...,Af, under the same assumptions as above, 
we derive from (3.21) 

g(k) -g n = a (k -in n )* N , n ® * N , n t N g(k) + g N . n u'[g(k) - g n \ + (i). 

Thus again, as in derivation of (4.15), we have by (4.14) that 

g(k) - g n = 0(k -iK n )$ n ® y N ,n^N 5(k) + (1), (4.29) 

where notation (4.16) was used. In order to find the expression for ^fN,n~^N g(k) we apply 
operation ^jv,n^jv to the above equality from the left, that gives 



*^ g W = l-^ + ° (1) 



(k- 

recalling (4.17) and the second equality in (4.22). Notice, that thanks to (3.13) and (4.25) 
the denominator here is different from zero, so we get by (4.29) 

g(k) = g n +a(k-iK n ) - *"®\ ; + o(l) (4.30) 

for k ~ in n for any n G [1, N]. 

4.3 Behavior of $(k) and ^(k) for k ~ m n 

Thanks to (4.30) and definitions (4.7) we can now derive the behavior of the Jost solutions 
in vicinities of the points of discontinuity. Thus, applying C N$jv(k) from the right to (4.30) 
we get, for k ~ in n 

$(x,k) = - gffi , +o(l), k~m„, n = l,...,Af, (4.31) 

1 - c n a(k-iK n ) 

recalling that ^^(k) is continuous at k = in n , and using equalities (4.18) and (4.21). 
Situation with ^(x, k) is a bit more complicated due to the pole behavior of ^jy(k) at these 
points. Thus, applying ^^(k)!? n from the left to (4.30) we get, recalling (4.19) and (4.22) 

We see that reductions of the relation (4.30) for the Green's functions describe the 
asymptotic behavior of the Jost solutions in vicinities of their discontinuities. Notice that 
this behavior is in agreement with (4.18) and (4.19) thanks to properties of the function 
a(k) in (3.12)-(3.14). Additional relations on the discrete values of the Jost solutions follow 
by means of the reduction of (4.28). Say, applying C n&n,u to this equality from the right 
we get by (4.16) 

l=m-\-l 
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where we have introduced the constants 

c mn = (^ m ^N^N,n), m,n = l,...,Af. (4.34) 
Definition (4.34) thanks to (4.1), (4.8), and (4.17) can be written in the different forms 

C mn = (^N,m^N Qm ^N$N,n) = (* m U'$N,n) , (4-35) 

where, say, the last equality is a short form for 

c mn = J dxu'{x)y m {x)® N ^{x). (4.36) 
We see that the diagonal terms coincide with (4.21): 

C nn = C n - (4.37) 

In contrast to the diagonal case, the off-diagonal constants looks to need more strict condi- 
tions on the perturbation u'(x). Indeed, as follows from (4.24) the product x F m (x)$ A , n (x) 
in (4.36) can grow exponentially, so u'(x) must decay at space infinity faster than some 
exponent linear in x, in order that integrals in this equality converge. 

Eq. (4.33) enables to establish relations between different discrete values, generalizing 
(2.30). Thus, summing it up with the matrix T> n j we cancel the term in the l.h.s. thanks to 
(2.30), getting 

N JV n+Af 

3 = 0. (4.38) 



E*»^+E E 

n—l n—ll—m+1 I 



The m-dcpcndencc of the second term is irrelevant: subtracting from (4.38) the same equa- 
tion with m->m-lwc get, due to (2.30), only a trivial consequence of (4.34) 

Af 

J2c mn T> nj =0. (4.39) 

n=l 

Let now m = J\f in (4.38). This gives a generalization of the first equality in (2.30) for the 
perturbed case: 

($ 1 (x),...,$ A ,(x))2? = 0, (4.40) 

where 

Af i n-1 



= V nj +T— E °nl V l3 = V nj ~ T^T" E ^ V lj = 

n l=n " 1=1 

s ih-( v « + £ CniVij )- (4 - 41) 



l=n+l 



the second and third equalities following from (4.39) and (4.37), correspondingly. Since for 
the zero perturbation, u'(x) = 0, the matrix V reduces to the matrix T>, that we choose to 
be TP, it is natural to assume that also the matrix V is TP. 
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5 Inverse problem 

5.1 <9- derivatives of the Green's function and Jost solution 

In order to find the departure from analyticity of the Green's function we notice first that 
by (3.12) in the sense of distributions 

= sgnky 
<9k 27rik ' 

Therefore, for the d-bar derivative of the Green's function £/(k) in a neighborhood of a point 
k <~ in n we have thanks to (4.30), 

dg(k) _ sgnk^ $ n ®*« 0(1) f52) 

<9k 27rz(k + iK„) [l-c n a(k-i Kn )} 2 

that is also valid in the sense of distributions, where the expression in the r.h.s. is locally 
integrable. Then the derivative of (4.4) sounds as 

ok 27r c*k 

where we used (3.10), (4.1), and (4.5). Thanks to (4.6) and the second equality in (4.7) this 
can be written as 

= $Jv( _k) ^ * ( -k) + M^'^- (5.3) 

In the r.h.s. of this equality we have a product of two distributions. Nevertheless, it is valid 
as thanks to (5.2), (2.28), (4.32), and (3.13) for the coefficients of the singular terms in both 
sides we have 

$n ® *„ = [1 - c„a(k -iKn)]$ Nin ® tf „ + [Q N<n +a(k -m™)^ ® *^„] ® * n) 

that is an identity, because for the first term in the second brackets we can use Q N, n U'$ n = 
®n — &N.n, that in its turn follows from (4.14) and (4.16), and in the second term we use 
the second equality in (4.22). Thus under the standard assumption of unique solvability of 
the integral equation (4.4) and applying C n$n{— k) to it from the right we derive 

^ = ^g ( k)^ ( -k)^ ( -k), 

that thanks to (4.6) and (4.7) gives finally 



dQ{k) _ sgnkft 
<9k " 2tt 



*(-k) ® tf(-k), (5.4) 



that due to the above consideration is valid in terms of distributions. 
From (4.7) and (5.4), recalling that $jv(k) is analytic, we get 



<9k 

with 



d*(x,-k) =$( ^ _k)r(k), (5.5) 



-(k) = ^(*(-E)1c JV * w (k)) > (5.6) 
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which, by using (4.7) and (4.6), can be rewritten as 

r(k) = !S^^ N( _^2 N g ik )t N ^ N (k)). (5.7) 

From (5.6), thanks to the behaviors of ^(k) derived in (4.32), and thanks to (4.34), (3.13), 
we obtain that r(k) is a singular function for all k = in n , n = 1, . . . ,Af and has the following 
singular structure: 

Therefore, thanks to the behavior of 3>(k), sec (4.31), we derive from (5.5), , that 

-= toejEB - + 0(1), k-i*,, (5.9) 

9k 27ri(k + ?k„)[1 — c n a(k — zk„)] 2 

i.e., this derivative has integrable singularities in all points k = in n . 
5.2 Integral equation of the Inverse problem 

Thanks to (5.9) $(x, k) is smooth enough and its 9-derivative does not need any regular- 
ization. But in order to derive the integral equation of the Inverse problem it is necessary 
to take the nonstandard normalization (4.11) into account. Thanks to (5.5) we have that 

$(x,k)e ika;i+k2a;2 = X(x,k) - - [ 4-^e lk ' a:i+k ' 2a:2 r(k')*(a;,-k 7 ), (5.10) 

7r J k — k 

where X(x, k) is an entire function of k. The integral term decays when k — > oo, so by 
(4.11) X(x, k) is a polynomial of order Nb with coefficient i Nb of the highest power. The 
new potential u(x) is defined then in analogy to (2.35) by 

u(x) = -2 lim (ik)- JVi > +1 a El (e ikxi+k2 X2 <$>(x,k)) . (5.11) 

k^oo 



So, as follows from (5.10), 



u(x) = -2 lim (ik)- Nb+1 X Xl (x,k). (5.12) 

k— t-oo 



Thus in order to get a closed system of equations of the Inverse problem we need to determine 
X(x, k) in terms of the Jost solutions itself. In particular, values X(x,in n ) must be expressed 
in terms of the values $„(x) of the Jost solution by means of the asymptotic relations (4.31). 
Such expression will allow to use condition (4.40). Because of the asymptotic behavior in 
(5.9) the limit k <~ in n of the integral term in (5.10) needs a special consideration. For this 
aim we use the following [18] definition of the principal value distribution: 

M 2 k'.f(k')sgnk^ , f <Pk> /(k')sgnl4 

p - v -y |k'p[i_ ca( k')] 2 ~iTo J |k'p[i- ca( k')] 2 ' {bAS) 

|k'|>e 

where /(k) is an arbitrary test- function. 

Statement. Let k, k' € C, a(k) be defined in (3.12) and c an arbitrary real constant. Then 
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1. the distribution (5.13) can be given equivalently by any of the next two equalities 

f d*k' /(kQsgnk^ 1 f 2 I , /(kQ-ZC-k 7 ) 

, /• d 2 k / [/(k / )-g(l-|k / |)/(0)]sgnk^ 
~y |k'|2[l-co(k')] 2 ' 1 ' 

2. in terms of distributions the following asymptotic is valid: 
sgnk^ _ sgnk^ 



k'(k'-k)[l-m(k')] 2 '|k'| 2 [l-ca(k')] 2 

Proof. Equality (5.14) results from (3.13) thanks to sgnk'j} and does not need any regular- 
ization. This equality proves that definition (5.13) is meaningful. Now (5.15) follows from 
(5.14) again thanks to (3.13) and sgn-function. In order to prove (5.16) we write 

/■ _rf 2 k'/(k')sgnk^ f d^k' [/(k')-fl(l-|k'|)/(0)]sgnk^ 

J k'(k'-k)[l-ca(k')] 2 J k'(k'-k)[l -ca(k')] 2 

+ /(0) / _ d2k '^ , 

J {<1 k'(k'-k)[l-ca(k')] 2 ' 

where the first integral now admits limit k — >• 0. In the second term like in derivation of 
(5.2) we assume that c is different from zero and use that thanks to (5.1) 

sgnk^ 2ttw 1 



k 7 [l-ca(k)] 2 ~ c k 'l-ca(k')' 

Then 

rf 2 k' /(k')sgnk^ _ fd^k' {f(k')-6(l-\k'\)f(0)}sgnk' 3 



« + 



k'(k'-k)[l-m(k')] 2 J |k'| 2 [l-ca(k')] 2 

+ /(0)— / rf 2 k'^ = l +0 (i), k ~0. 

W c 7 k'-kak'l-ca(k') V h 



k' |<i 



Thanks to (5.15) the first term is the principal value distribution (5.13). For the second 
term thanks to the Cauchy-Green formula for | k | < 1 we have 

f 2 , 1 <9 1 _ 1 f , 1 7T 

7 k' - k ak 7 1 - co(k') ~ 2^ / (k'-k)[l-ca(k')] l-co(k)" 

k'|<l I k' 1=1 

Now we can perform the limit in the contour integral, getting 

r , 1 r dk' _ f da 

kTo y dk (k'-k)[i- ca (kO] " y k'[i- ca (k')] j i - ca/ir - 

I k' 1 = 1 I k' 1 = 1 
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-2m 



log(l - c), 



where (3.14) was used. Summarizing, we have that 



/ 



d 2 k' /(k')sgnk(ft 
k 7 (k'-k)[l-ca(k')] 2 



p.v. 



d 2 k' /(k')sgnk^ 
|k'| 2 [l-ca(k')] 2 



-/(o) 



27T 2 i 



log(l - c) 



1 - ca(k) 



o(l), k~0. 



that proves (5.16). The case c = follows then by the limiting procedure. ■ 

Now we can get values of X(x, k) at points k = in n , n £ [1, N], in terms of the discrete 
values (4.16) of the Jost solution. Thanks to (5.9) the singular behavior of the integrand in 
(5.10) at k ~ in n is exactly like in (5.16), once shifted k and k' for —in n . We define 

,2 k / 



r d 2 k' 

P J k' -in. 



.ik'xi+k' X2„(lrl 



r(k')$(x,-k') 



lim 

e->0 



d 2 k! 
k' 



r(k')$(x,-k'). 



(5-17) 



I k' —%K n I >e 



so that by (5.10) and (5.16) 

<S>(x,k)e lkxi+k2x2 =X(x,in n )- 

-Ip.v. / f k e ik ' Xl+k " x *r(\*!)<f>(x,-k!) 
it J k —in n 



+ <f> n (x)e 



-K n (x) 



log(l - cj 



1 - c n a(k-i/c„) 



+ o(l), k~i« n , 
(5.18) 

where notation (2.4) was used. Now by (4.31) we get that the undetermined terms on both 
sides cancel out and we arrive to 



X(x,iK n ) = <p n (x)e K "^ + 



+ -p.v. 

7T 



d 2 k' 



ik'rri+k' x-2 (i.! 



where we introduced 



<Pn{x) 



k' — in. 
log(l - cj 



r(k')$(x,-k'), 



$„(z), ,n = l,...,Af. 



(5.19) 



(5.20) 



log ( 1 — ■ Cji ) 

Notice that the ratio — is positive thanks to (4.25). Let in analogy to (2.28) 

(p(x) = {(p 1 (x),...,(pjs(x)}, (5.21) 
denotes a AA-row of rescaled solutions of the heat equation. Then (4.40) is written as 

ip(x)V = 0, (5.22) 

where (see (4.41)) 

-c„ 



Vnj (l-cJlog(l-cJ 



l=n+l 



V n + c ni v ij), n = l,...,AT, j = l,...,N b , (5.23) 



that gives Nb equations on Af = N a + Nf, unknowns ip n (x) 
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5.3 Closure of the Inverse problem 

Assuming that the values of the polynomial X(x, k) at points k known, it can be 

written as 

X(x^) = J2x(x,i Kn )f[ *i±i* (5.24) 

n=l m=l m " 

that, apparently, is a polynomial of power J\f — 1, in contradiction with the properties of this 
polynomial mentioned after (5.10). In fact, in order the get an expression for X{x, k), which 
explicitly displays the required properties, we need to take into account that the values 
X(x,in n ) are not independent. In order to formulate these conditions, we denote by 

Sj = Sj («!,..., K/v) = E K m--- K no S = 1) 

l<n 1 <...<n i <7V 

the symmetric polynomials of their arguments and by sj n ^ — • • • , k^, ■ ■ • So™" 1 = 1 

the symmetric polynomials with omitted variable K n . Then sj™^ = X^j=o( — R n)^ s i-j an d 

(«m + ik) = ^ (ik^- 1 - 1 ^ = ^ (ik)^" 1 -' (5.25) 
m=l z=o z=o j 

Thus, (5.24) thanks to (2.9) takes the form 

X(x, k) = (-1)" 1 £ (i k)" ^ £ A(x, i« n )(-« n y 7n . (5.26) 

Z=0 i=0 n=l 

In order to cancel extra powers all terms with I = 0, . . . , N a — 2 must be zero, where by (2.2) 
N a = M —Nf,. Choosing first / = we get that the term with j = in the sum over j equals 
zero. Then both sums by I and j goes from 1. Now, in order to cancel the term with I = 1 
we have to put that term of the j-sum with j = 1 equal to zero. Repeating this procedure 
up to I — N a — 2 we derive that all terms up to j = N a — 2 equal zero. Coefficient of the 
term with I — N a — 1 must be equal to i Nb in order to obey the normalization condition for 
polynomial X(x,k) (see remark after (5.10)). All together this gives that the values of this 
polynomial must obey N a relations 

M 

Y,X(x,iK n W„'Yn = (-l) Nb 6i,N a -i, l = 0,...,N a -l, (5.27) 
and representation (5.26) takes the form 

N b 

X(x,k) = J2(iV Nb - l Si+ 

1=0 

N b I AT 

+ (-V*- 1 £(* E s i-i E X ( x > ^n)(-Kj j+JV »- 1 7„. (5.28) 

1=1 j=l n=l 

Relations (5.27) due to (5.19) give 
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= (-1)^,^-1, l = 0,...,N a -l. 
Thanks to (2.9) it is easy to check that 

e ^zt - -n-*f-\i n ^ip. » = o, . . . , jv 0> 

n=l n n=l ™ 

so we can write relations (5.29) in the form 

<p(x) ie - K V = ( 0, . . . , , (-!)"» ) + v(x), 

JV„ 

where we introduced iV a -row 

v(x) = {v^x),. . . ,v Na (x)}, 

f , 2l , ^k')'- 1 e ik ' a:i + k ' 2a;2 r(k')$(a;,-k 7 ) 
vi(x) = - — 7z p.v. / eTk 



Z7T 



(5.29) 
(5.30) 
(5.31) 

(5.32) 
(5.33) 



for Z = 1, . . . , -/V a + 1, being the value Z = Af a + 1 needed below. Together with (5.22) these 
equations give just N equations to determine N discrete values of the Jost solutions, that 
close the formulation of the inverse problem. In order to reconstruct the perturbed potential 
we notice that by (5.12) 



u{x) = -2 lim (ik)- Nb+1 X Xl (x,k), 

k— >oo 



so that since (5.28) we get 



A/" 



U{X) = 2(-lf*d Xl J2 < a lnX{x, iK n ), 



n=l 



that in its turn can be written thanks to (5.19) and (5.30) as 

u(x) = 2(-l)^a ci (J2 C°7„e-^ (x V„(z) - v Na+1 (x)) , (5.34) 

where v^ a+ i(x) is defined in (5.33). 

Finally, we have to exclude X(x,k) from the equation of the Inverse problem (5.10). 
Inserting (5.19) in (5.24) we get 

X(x,k) = J2Mx)e- K " {x) n + 

n—1 m.—l 



l 



—p.v. /(ike 

7T 



2 lr' J- k ' ^l+k' x 2 



1 T r ^rn ^ k 



Kk')^-k')El7Z— II 



i k %K n K m K n 

n=l m=l 
ra^n 



From the trivial identity 



1 -Q Atm + »k _ 1 

i k Kn k k 

n—1 m— 1 



j-j /t„ + i k _^ 

,n=l K " + 1 k 
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we have that (5.10) is reduced to the form 

n=l m=l Km K " 

+ W CT eik ' Ki+k ' 2Kv(k ' )$( ^- F) n^^- (^35) 

Thus the system of equations of the inverse problem is given by (5.35), (5.22), and (5.31). 
In terms of solution of this system the potential is reconstructed by means of (5.34). 



6 Time evolution 

In [17] we proved that the multisoliton potential of the heat operator constructed by means of 
the Darboux transformation coincides with the potential derived by the r-function approach. 
Thus we can use the standard scheme of this approach, see [25], where the time evolution 
of the pure TV soliton solution of KPII is obtained by adding a time term to K n (x) defined 
in (2.4), i.e., by writing 

K n (x, t) = K n X\ + K n X2 — An n t. (6.1) 

Under this substitution Un(x) given in (2.10) obtaines time dependence and becomes the N 
soliton solution of the KPII equation (1.1). Correspondingly, the time evolution of the Jost 
solutions (2.12) is given then by the equations 

d t $ N (x, k) = -a£$ n (x, k) + 4« k 3 $jv(x, k), (6.2) 
d t ^ N (x, k) = * N (x, k)X< - Ai k 3 *at(x, k), (6.3) 

where An{x, 9) is given in (1.3) with u substituted by un and where we used notation (2.40) 
for the action of operator A(x, d) and its dual. Time evolution (6.1) can be reformulated 
equivalently preserving definition of K n (x), but including time dependence in the spectral 
data of the N soliton solution: 

(64 » 

^^ = -4 K lV nm +(Va) nmi n = l,...Af, m = l,...,N b , (6.5) 

where a is an arbitrary (possibly time depending) Nf, x matrix corresponding to the re- 
mark that the iV-soliton solution is fixed not by this matrix, but by a point on a Grassmanian. 
If we choose the matrix V in some way, then the matrix a must be fixed correspondingly. 
Say, if we choose (without loss of generality) that V nm = <5„, m for all m, n = 1, . . . , Nf, then 
a = diag{4K?,...,4K 3 v J. 

Let us notice that above results for the N soliton solution are independent of the choice of 
the sign in the bottom limit of the integral term in (1.3), as thanks to (2.10) and asymptotic 

behavior (2.25) of the r-function, we have that J dx\ d% 2 UN(xi, x 2 ) = 0. This situation 

— oo 

is changed when we consider the perturbed solution u(x). Understanding time dependence 
in all terms of (1.8) and taking that both u(x) and un(x) obey KPII equation (1.1) into 
account, we get that perturbation u' must obey 

(u' t - 6u'u' Xl - 6{u'u N ) Xl + u' XlXlXl ) Xl = Su' X2X2 . (6.6) 
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Again, using that Un(x) exponcntialy decays when x\ — > ±00, we see that asymptotically 
v! is governed by the KPII equation, more exactly by its linear part. Thus according to 
an analysis analogous to that performed in [26] for the KPI equation, the evolution form of 
(6.6) is given by 

u 't ~ 6u'u' Xl - 6(u'u N ) Xl + u' XiXiXi = -3 J dx' lU ' X2X2 {x' 1 ,x 2 ), (6.7) 

too 

where too denotes sign of the infinity limit. Thanks to the above discussion, the perturbed 
solution u(x) must satisfy the evolution equation 

/■Xl 

ut(x) - 6u(x)u Xl (x) +u XlXlXl (x) = -3 / dx' 1 u X2X2 (x' 1 ,x 2 ), (6.8) 

J too 

so that for generic intial data the time derivative of the solution is discontinuouse at t = 0. In 
order to get this form of the evolution equation by means of the compatibility condition (1.4) 
we have to modify, correspondingly, the definition of operator A in (1.3) correspondingly: 

Xl 

A(x) — Ad Xl — 6ud Xl — 3u Xl — 3 J dx[ u X2 (x' 1 ,x 2 ). (6.9) 

too 

In order to find out the time evolution of the Jost solution $(x, k) we notice that thanks 
to (2.28), (2.29) and (6.2), (6.3) 

d t $N,n( x ) = -A^$N,n( x ) + 4k^$at,„(x), (6.10) 

d t ^N,n( x ) = ^nA^-An - 4/4*iv,n(a;), (6.11) 
n = 1, . . . ,J\f, so that by (3.1) 

dtQ N = [QN,A N }. (6.12) 
Now differentiating (1.9) by t we use that (6.7) can be written thanks to (1.4) and (1.8) as 

u' t = [£ N ,A N ]-[£,A], (6.13) 

that by means of (6.2) and (6.12) enables derivation of the standard relations 

dt®(x, k) = -~%<b{x, k) + Ai k 3 <!>(x, k), d t ^{x, k) = V(x, k)^4 - Ai k 3 *(a;, k). (6.14) 

Then derivative of (5.5) by t gives 

<9 t r(k, t) = Ai(k 3 +k 3 )r(k, t), (6.15) 

so that from the singular behavior of r(k), $(x,k) and ^(x, k) (see (5.8) and (4.31), (4.32)) 
we derive 

fe=0, (6.16) 
a t $„(a;) = -^$„(x)+4 K 3 $ n (x), for n = 1, . . . , W, (6.17) 
&t* n {x) - * n (xj% - 44*nW- (6-18) 
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Thus by (5.20) 

dt<p n (x) = -~X<p n (x) +4K% l <p n (x), for n = 1, . . . , Af, (6.19) 

and following the same procedure as in derivation of (6.14) we get by (6.13), (6.10) and 
(6.18) for the time derivative of constants c mn defined in (4.35) 

dc 

= 4(nl - K 3 m )c nm , n,m = 1, . . . ,Af . (6.20) 
Now from (5.23), thanks to (6.5), (6.16) and (6.20), we derive 
(JbD ' — ' — 

—^=-4K 3 n T> nm + (T>a) nm , n=l,...Af, m = l,...,N b . (6.21) 



7 Conclusion 

We developed here a modification of the 1ST that enables to consider solutions of KPII of the 
kind (1.8). The study of a heat operator (1.2) with such potentials, having ray behavior on 
the x-plane, cannot be considered, strictly speaking, a "scattering problem" . Nevertheless, 
in terms of this modification we succeeded in introducing Jost solutions and scattering data 
and we proved that the mentioned asymptotic behavior of the potential results in specific 
discontinuities and singularities of these objects with respect to the spectral parameter 
k. Consequently, we gave a corresponding modification of the Inverse problem and time 
evolution of the scattering data. The proof of the unique solvability of the Direct (1.9) and 
Inverse problem (5.22), (5.31), and (5.35) was not our aim here. The whole construction is 
based on the assumption of the unique solvability of (1.9) for the class of potentials under 
consideration. While the corresponding proof must be a rather straight generalization of 
the case of a decaying potential (see [4]- [8]), since we have boundedness property (1.11), 
the investigation of the Inverse problem is much more involved and may need more details 
on properties of the scattering data. Here we mention only that (5.22) and (5.31) give a 
system of linear equations on the Af-row <p(x) in (5.21). Relation (5.22) can be considered a 
generalization to the perturbed case of relation (2.30) of the pure soliton case. If, in analogy 

to (2.6), there exists a matrix V (with at least one nonzero maximal minor) such that 

T>'T> = 0, then (5.22) means that there exists an N a row <p(x) such that 

<p(x)=<p{x)V'. (7.1) 
Then, for the unique solvability of relation (5.31) it is necessary that the determinant 

t'(x) =detT>'je- K V' (7.2) 

is different from zero. Relation (7.2) is analogous to the second representation for the pure 
soliton r- function in (2.11). Correspondingly, we inherit from the pure soliton case the 

condition that the matrix T> is TNN, or TP. Taking (4.25) into account, we get by (5.23) 
that the matrix T> n j + Y^iLn+i c nl ^ih where the matrix c„; is constrained by (4.39), must 
have the same property. It is also necessary to mention that the diagonal elements of the 
matrix c are related to the singular behavior of the spectral data r(k), see (4.37) and (5.8). 
Say, if r(k) = we have that all diagonal elements are zero. In this case the Inverse problem 
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(5.22), (5.31), and (5.35) supply us with another N soliton solution. This solution is given by 
a r-function (7.2) parameterized by the matrix T> n j + Y^iL n +i c «' ^ij instead of the matrix 
T> n j. A formulation of conditions for this matrix to be TP, or at least TNN, as well for the 
above mentioned unique solvability, in the generic case, of the Direct and Inverse problems 
are open at this moment. 
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